Metastable criticality and the super Tonks-Girardeau gas 
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We consider a ID Bose gas with attractive interactions in an out-of-equilibrium highly excited 
state containing no bound states. We show that relaxation processes in the gas are suppressed, 
making the system metastable on long timescales. We compute dynamical correlation functions, 
revealing the structure of excitations, an enhancement of umklapp correlations and new branches 
due to intermediate bound states. These features give a clear indication of the attractive regime and 
can be probed experimentally. We observe that, despite its out-of-equilibrium nature, the system 
displays critical behaviour: correlation functions are characterised by asymptotic power-law decay 
described by the Luttinger liquid framework. 



Quantum many-body systems exhibiting strong corre- 
lations are of prime interest from both theoretical and ex- 
perimental perspectives as they lead to the breakdown of 
the simple single-particle excitation picture. Indeed, col- 
lective modes reign over the low-energy sector, signalling 
the presence of quantum critical behaviour [TJ and refer- 
ences therein]. An example of such a system is the ID 
Bose gas [2] which due to recent experimental progress 
can be realised and probed in and out of equilibrium [3HS] 
with tuneable interaction strength [7]. 

An interesting question is: how strong can the cor- 
relation become? Increasing local repulsive interactions 
ultimately leads to the Tonks-Girardeau gas (TG) [EJ|9], 
which corresponds to infinite interactions and seems to 
set the limit of correlation strength in the ID Bose gas. 
However even stronger correlations are possible if we con- 
sider a specific metastable state of a gas with attrac- 
tive interactions, the super Tonks-Girardeau gas (sTG) 

Attractive point-like interactions drastically modify 
the dynamics of the gas as they allow the formation 
of bound states. In fact, the ground state is a macro- 
scopically large molecule [16] , whose dynamical responses 
can be computed exactly [T7] . Instead we consider here 
a completely out-of-equilibrium state of the system in 
which there are no bound states and all particles are in 
a low- lying gaseous state. Contrary to the ground state 
this state is characterised by a finite energy per particle 
in the thermodynamic limit. The structure of excitations 
around it is similar to that found in fermionic systems 
with the addition of bound states. 

Relaxation in this system thus goes hand-in-hand with 
the formation of bound states. However, as we will see, 
local fluctuations are not efficient in bringing particles to- 
gether, greatly suppressing the rate of formation of bound 
states. This makes the system sufficiently long-lived to 
be experimentally observable. 

Our starting point is the Lieb-Liniger Hamiltonian 
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where N is the number of particles and c parametrizes 
interactions. Throughout the Letter we focus mostly on 
the attractive regime (c < 0). We address the correla- 
tions explicitly by studying the density-density correla- 
tion function defined by its Lehmann representation 

S(k,w) = f E \(MPk\sTG)\ 2 5 (to — E\ + E sTQ ) , (2) 
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where A labels states in the Hilbert space pk is 
the Fourier transform of the density operator {pk = 
I Hq A) and is the energy of a state |A). |sTG) 
is the sTG state and its explicit definition is provided af- 
ter we discuss the eigenstates. The density-density cor- 
relation can be measured using Bragg spectroscopy [T8] . 

The sTG gas can be experimentally realised by em- 
ploying a confinement induced resonance [SJ [7] . Quench- 
ing the magnetic field over the resonance changes the 
sign of the interactions between particles. Starting in the 
ground state of the TG gas, the interaction parameter is 
quenched to the other side of the resonance. Quenching 
to c = — oo exclusively populates the sTG state. Quench- 
ing to a finite but large negative value of c leads to a 
combination of states with the sTG state being predomi- 
nant, all other states having an amplitude suppressed by 
at least a factor |c| -2 . 

The Letter is organised as follows. We start with the 
description of the eigenstates, first generally for attrac- 
tive interactions and later specifically around the sTG 
state. We then discuss the results for the density-density 
correlation in momentum space and real space, display- 
ing their metastable quantum critical features. 

Eigenstates. The Lieb-Liniger Hamiltonian (fil) ex- 
hibits an exact solution by the Bethe Ansatz [2]7^ The 
N-particle wave function is given by a linear combination 
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of plane waves 
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with Ap = (-1) [P] e^>»' gn{xi - Xh) ^- Xp ^ and the 
two-particle phase shift </> (A) = 2 arctan (A/c). By im- 
posing periodic boundary conditions rapidities {Aj}^ 
get quantised and fulfil a set of coupled non-linear equa- 
tions (Bethe equations) 
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The quantum numbers {Ij}^ =1 are integers (N odd) or 
half-odd integers (N even) and completely specify the 
eigenstate. The Hilbert space of Hll is spanned by 
different choices of sets {Ij}. The momentum of the 
eigenstate is P\ = Ylj=i Aj , whereas the energy equals 

The structure of the Hilbert space is different depend- 
ing on the sign of interactions. Characteristic for repul- 
sive interactions is a Pauli-like principle: quantum num- 
bers must be mutually distinct and all solutions to the 
Bethe equations are real. Attractive interactions drasti- 
cally change this by allowing complex solutions [16] . Ra- 
pidities then form regular patterns in the complex plane 
(strings) with exponentially small corrections in the sys- 
tem size which for the system size considered here are 
completely negligible. They are symmetric around the 
real axis and their imaginary parts are separated by \c\. 
These complex states are naturally understood as bound 
states and should be viewed as independent, stable par- 
ticles. The non-zero imaginary part of rapidities causes 
exponential decay of the wave function for an increasing 
separation of two particles within a string with a charac- 
teristic distance |c| _1 [TB] . 

The ground state is formed by a single bound state 
(TV-string). All other eigenstates can be constructed by 
two basic operations, either by breaking them and/or by 
giving them momentum. In the attractive regime, the 
Pauli principle holds between strings of the same length. 
The state classification is then straightforward [T7] . 

sTG gas. For a system with N particles we explicitly 
define the sTG gas by the following choice of quantum 
numbers 

I fTG = _N+l +jt j = 1; ... ;iV . (5) 

The quantum numbers are exactly the same as for the 
ground state of the repulsive gas of N particles and form 
a Fermi sea (Fig. [I]) . 

In order to compute correlations we first need to clas- 
sify the excitations around the sTG gas. One class of 
excitations contains the particle-hole excitations known 
from the repulsive gas [2] . We can again distinguish type 
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FIG. 1. Pictorial representation of quantum numbers for sys- 
tem with N — 8 particles and in a) the sTG state, b) an 
excited state with one particle-hole excitation c) an excited 
state with one 2-string. 



I and type II modes which can be combined to create 
multi particle- hole excited states (see Fig. [TJd). These 
particle-hole excitations lead to a modified version of cor- 
relations as compared to the repulsive gas. Attractive 
interactions induce smooth changes to the particle-hole 
contribution to the correlation as one tunes the variable 
1/c through zero. 

Other classes of excitations involve bound states. They 
are directly understood as processes of binding single 
particles in molecules (see Fig. [l^) and are ultimately 
responsible for the instability of the gas. These classes 
lead to new branches of correlations which we discuss in 
the next section. Before that however we examine the 
simplest (and most probable) process: the formation of a 
2-string (bound state with 2 particles) in order to quan- 
tify the timescale for stability of the gas. 

Consider adding a localised impurity potential to the 
otherwise isolated system. This perturbation connects 
the sTG state with all the states at the same energy. The 
timescale can be directly computed using Fermi's golden 
rule. Assuming that the external potential is weak (we 
set the magnitude of the potential to be V = O.lep ) 
for a gas of 133 Cs atoms of ID density n = 10 6 m _1 with 
c/n -8 we get r 2 _ str 2hTiV 2 S(uj 0")/2m ~ 1 s _1 . 
S(u = 0~) is the density-density correlator summed over 
all momenta (see Eq. |8| and taken at the slightly neg- 
ative value of the energy so that only bound states con- 
tribute. For larger attractive interactions and weaker 
perturbations the timescale is much larger. The typi- 
cal timescale of experiments with cold atomic gases are 
shorter (~ 1(T 3 - 10" 1 s) than I^ 1 tr and thus the sTG 
gas is experimentally stable. The rate of formation of 
bound states is approximately only one order of magni- 
tude larger than the rate of 3-body collisions [19]. This 
shows that the sTG gas can be viewed in practice as be- 
ing as stable as the repulsive gas. 

A physical picture of this statement on stability is as 
follows. For a 2-string to form, two particles must be a 
distance |c| _1 apart. This is however highly unlikely since 
the initial non-local pair correlation function (S(x)) ex- 
hibits fermionic behaviour: at small distances its value 
is small. This shows that possible decay channels are 
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TABLE I. Levels of saturation of the f-sum rule. All com- 
putations were performed at unit filling (N/L — 1) and for 
N = 128 particles. 
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not really open (for c = — oo they are in fact completely 
closed since the rates from Fermi's golden rule vanish). 
We further analyse the behaviour of S(x) in the next sec- 
tion. Further evidence for the stability of the sTG gas 
is provided by numerical computation of the compress- 
ibility of the gas which is positive for large attractive 
interactions [TO] . 

Density- density correlation function. We directly 
compute the Fourier transform of the density-density cor- 
relation function (Eq. [5| by adapting the method used 
previously for the repulsive gas [20 . The matrix elements 
of the density operator (form factors) are known exactly 
through the methods of Algebraic Bethe Ansatz [21 j. Eq. 
([2| demonstrates the computational method. The ABA- 
CUS algorithm [22] was used to recursively explore the 
Hilbert space looking for states that contribute the most 
to the correlation. Upon summing contributions we ob- 
tained highly accurate results for to the density-density 
correlation function for N = 128. The exact f-sum rule 
identity 

j°°^S{k^ = ^k\ (6) 

gave a quantitative check of the calculations (Tab. 

Momentum space. The full dynamical correlation is 
plotted in Fig. [2] As the sTG state takes the form of 
a Fermi sea, its low-lying type I and II excitations have 
a linear spectrum and this sector of the sTG gas falls 
into the Luttinger liquid universality class, with Lut- 
tinger parameter 0.5 < K < 1. The behaviour of the 
correlation along the edges of support of a single particle- 
hole excitation agrees in fact with the predictions of non- 
linear Luttinger theory [23]. For values of K > 1 (re- 
pulsive interactions) there is a singularity along the type 
I mode (upper threshold) and smooth vanishing of cor- 
relation along the type II mode (lower threshold). For 
0.5 < K < 1, inversely to the repulsive case, the correla- 
tion is smooth along the particle mode and diverges along 
the hole mode. This shift of correlation weight towards 
the lower threshold is exactly what we observe (see Fig. 
[2]and[3|. 

Two other interesting quantities namely the static cor- 
relator and the dynamical autocorrelator 

S(k) = J? 00 S(k,u>)%, (7) 
S(u) = iE k S(k,u), (8) 




FIG. 2. (color online) The density-density correlation in mo- 
mentum and energy space for c — —8, c — —16 and c— —64. 
The negative energy parts of the plots were rescaled by factors 
10, 100 and 2.5 x 10 5 respectively to make the string contri- 
bution easier to see. The discontinuity in the c = — 8 plot 
around uj — is an artefact of this rescaling. As the interac- 
tions become more attractive the correlation weight spreads 
uniformly between the lower and upper thresholds of the sin- 
gle particle- hole continuum, just as for the TG gas. At the 
same time the contribution from bound states is suppressed 
in value and moves to lower energy. 
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FIG. 3. (color online) Fixed momentum cuts through the 
dynamical structure factor. The shift of the correlation weight 
towards the lower threshold and the contribution from the 
bound states for uj < are clearly visible. For comparison 
the c = 16 ground state correlation, where the singularity is 
at the upper threshold, was also plotted. 



are plotted in Fig. [4] The shift of the correlation weight 
towards the lower threshold leaves a characteristic trade- 
mark in the static correlator (Eq. [7| and Fig. [4| [TU] . 
For repulsive interactions, the static correlator around 
k = 2kp smoothly approaches its asymptotic value from 
below (S(k) — »• 1). Here at k = 2kp we observe a di- 
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FIG. 4. (color online) Top: Static structure factor. Attrac- 
tive interactions lead to a singularity at k = 2k f- Bottom: 
Dynamical auto correlator. For smaller attractive interactions 
an extended plateau develops for uj < 0. This clearly indi- 
cates the contribution to the density-density correlation from 
the bound states. 




x[L] 

FIG. 5. (color online) Distance dependence of the density- 
density correlation. Despite varying the interaction param- 
eter system is consistently fermionic in behaviour. Smaller 
attractive interactions increase Friedel oscillations. Inset: 
Correlations immediately approach the asymptotic behaviour 
predicted by Luttinger liquid theory. All curves converge to 



vergence of the correlator and a power-law tail above the 
asymptote. 

Moreover, for less attractive interactions, the dynam- 
ical autocorrelator S(cj) (Eq. [8] and Fig. [4| develops a 
plateau on the negative side of uj. This plateau is a clear 
signature of the attractive gas and of the existence of 
bound states. Together with the divergence at the lower 
threshold, these are the two smoking guns to look for 
experimentally. 

Real space. We now move on to real space and inspect 
the Fourier transform of the static correlator 



S(x) 



-ikx 



L 



S(k). 



(9) 



From the behaviour of S(x) for small x we can infer the 
effective statistics of the particles (Fig. [5|. We see the 
fermionic-like behaviour which is robust to changes in 
the interaction [15]. This is similar to the fermionization 
process in the repulsive ID Bose gas [9]. 

At large distances, results can be compared with Lut- 
tinger liquid theory which predicts [l] [24] the large dis- 
tance asymptote of correlations as a series 
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This formula is valid in the absence of bound states. 
These can be captured by treating them as mobile im- 
purities in the gas. The bosonization of the system then 
leads to a theory similar to that of the ferromagnetic Lut- 
tinger liquid [25] . However the existence of bound states 
does not have a strong effect on the spatial correlator giv- 
ing only subleading corrections which we will investigate 
in future publications. 



The pref actors A m in Eq. (10) are not universal. They 



depend on the microscopic theory and are connected with 
the scaling limit of form factors [26] . This in turn allows 
for the exact computation of prefactors in integrable the- 
ories where form factors are explicitly known. Fig. [5] de- 
picts the comparison between the correlations computed 
for N = 128 and the asymptotics (including prefactors) 
from [26] . The accurate match is yet another illustration 
of the applicability of the theory of Luttinger liquids in 
this out-of-equilibrium context. 

Conclusions and outlook. In this paper we considered 
the met ast able super Tonks- Girardeau gas and its corre- 
lations. We showed that attractive interactions modify 
the dynamical responses of the system in a clear, exper- 
imentally observable way. The system exhibits stronger 
collective behaviour with the majority of the density- 
density correlation carried by the type II mode. At- 
tractive interactions strengthen the umklapp excitations 
leading to an extended region of high correlation around 
2kp, which ultimately causes the divergence of the static 
correlator at k = 2kp- On top of this there are bound 
states with an extended region of correlation for uj < 0. 
These features provide experimentally clear signatures of 
the super Tonks- Girardeau gas, which should be experi- 
mentally accessible with current methods. 

We also showed that, despite attractive interactions 
and metastability, the super Tonks- Girardeau gas still 
displays the standard features of a quantum critical liq- 
uid. The system has a sector of excitations which falls 
into the Luttinger liquid universality class and the lead- 
ing long-distance asymptotes of correlations agree with 
the predictions of Luttinger liquid theory. By treating 
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bound states as impurities, a theory similar to the ferro- 
magnetic Luttinger liquid can be developed, allowing to 
study the dynamics of bound states in more detail. More- 
over, initial metastable quantum critical states can exist 
in other systems (e.g. multicomponent bosons, fermions). 
We will investigate these issues in future publications. 

We acknowledge useful discussions with H.-C. Nagerl 
and V. Gritsev. We thank SARA Computing and Net- 
working Services for access to the Lisa Compute Clus- 
ter. We gratefully acknowledge support from the Foun- 
dation for Fundamental Research on Matter (FOM) and 
from the Netherlands Organisation for Scientific Re- 
search (NWO). 



m.k.panfil@uva.nl 

[1] T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, 2004)_. 

[2] E. H. Lieb and W. Liniger,|Phys. Rev. 130, 1605 (1963)} 
E. H. Lieb, Phys. Rev. 130, 1616 (1963)| 

[3] T. Kinoshita, T. Wenger, and D. S. Weiss, Nature 440, 
900 (2006). 

[4] S. Hofferberth, I. Lesanovsky, B. Fischer, T. Schumm, 

and J. Schmiedmayer, Nature 449, 324 (2007). 
[5] E. Haller, M. Gustavsson, M. J. Mark, J. G. Danzl, 
R. Hart, G. Pupillo, and H.-C. Nagerl, Science 325, 
| 1224 (2009)1 
[6J S. Trotzky, Y.-A. Chen, A. Flesch, I. P McCulloch 
U. Sc hollwoeck, J. Eisert, and I. Bloch, |Nat Phys 8,| 
| 325 (2012)| 

[7] a 



[8] L. Tonks, |Pfys~Rev . 50, 955 (1936)| 
[9] M. GirardeauT p. Math. Phys. 1, 516 (I960)] 
[10] G. E. Astrakharch ik, J. Boronat, J. Ca sulleras, and 

S. Giorgini, [Phys. Rev. Lett. 95, 190407 (2005)| 
[11] M. T. Batchelor, M. Bortz, X. W. Guan, and N. Oelkers, 

J. Stat. Mech.: Th. Exp. 2005, L10001 (2005). 
[12] S. Chen, X.-W. Guan, X. Yin, L. Gua n, and M. T. 

Batchelor, |Phys. Rev. A 81, 031608 (2010) 
[13] M. D. Girardeau a nd G. E. Astrakharchik~ |Phys. Rev. A 



[14 

M 
[17; 



[is; 

[19 



[20 



[22 
[23' 



81, 061601 (2010) 



D. Muth and M. Fleischhauer, Phys. Rev. Lett. 105, 
150403 (2010). 
M. Kormos, G. Mussardo , 
Rev. A 83, 013617 (2011)] 
J. B. McGuire, J. Math. Phys. 5, 622 (1964). 



and A. Trombettoni, |Phys.| 



P. Ca labrese and J.-S. Caux, Phys. Rev. Lett. 98, 1 50403 
(2007): J. Stat. Mech.: Th. Exp. 2007, P08032 (2007). 
D. Cle ment, N. Fabbri, L. Fallani, C. Fo rt, and M. In- 

guscio, |Phys. Rev. Lett. 102, 155301 (2009)] 

D. M. Gangardt an d G. V. Shlyapnikov, |Phys. Rev. Lett.| 
90, 010401 (2003)| 



J.-S. Caux and P. Calabrese, |Phys. Rev. A 74, 031605 



(2006) 

N. A. Slav nov, Th eor. Math . Phys. 8 2, 273 (1990). 



J.-S. Caux, J. Math. Phys. 50, 095214 (2009 



3 (19 

mi 



A. Imambekov and L. I. Glazman, Phys. Rev. Lett. 100, 
206805 ( 2008) ; Science 32 3, 228 (2009). 
"[24] F. D M. Haldane, |Phys. Rev. Lett. 47, 1840 (1981)| 
[25] M. B. Zvonarev, V. V. Cheianov, and T. Giamarchi, 

|Phys. Rev. Lett. 99, 240404 (2007)] 
[26] A. Shashi, L. I. Glazman, J.-S. Caux, and A. Imam- 
bekov, Phys. Rev. B 84, 045408 (2011)} A. Shashi 
M. Panfil, J.-S. Caux, and A. Imambekov, Ph ys. Rev. B| 



[7] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998) 



85, 155136 (2012) 



